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bstract 

This work is devoted to analysis of the nonhnear dynamics of interrelated thermal and electro- 
magnetic perturbations in a superconductor caused by dissipative effects involved in the viscous 
motion of the magnetic flux. The structure and evolution of the stationary thermomagnetic waves 
propagating in the superconductor is studied based on a system of nonlinear equations describing 
dynamics of the instability development in the system studied. The apperance of these structures 
is qualititively described and the wave propagation velocity is estimated. The affect a thermoacti- 
vated flux creep to the structures of the nonlinear thermomagnetic wave propagation is studied. It is 
shown that taking into account a nonlinear dependency of the current density j on the electric field 
strength E would not qualitatively change the main results, since the character of the equilibrium 
state on the phase plane remains esentially the same. 

At present superconducting systems with high critical field strengs and current densites 
are widely implemented in various advenced technologies. Hovever, successful operation of 
the superconducting materials is only possible provided that special measures are taken to 
prevent a system from the thermal or magnetic breakage of superconductivity and transition 
to a normal (resistive) state [1]. For this reason, one of the main problems in the investigation 
of properties of superconductors is that of predicting the superconducting state breacage 
caused by dissipative and nonlinear effects related to viscous motions of the magnetic fiux. 
This explains a considerable interest in the study of dissipative and nonlinear effects in 
superconductors that has apose in recent years. 

This work is devoted to analysis of the nonlinear dynamics of interrelated thermal and 
electromagnetic perturbations in a superconductor caused by dissipative effects involved in 
the viscous motion of the magnetic fiux. The structure and evolution of the stationary 
thermomagnetic waves propagating in the superconductor is studied based on a system of 
nonlinear equations describing dynamics of the instability development in the system is 
studied. 

The evolution of thermal and electromagnetic perturbations in a superconductor is de- 
scribed by a set involving the thermal conductivity equation [2-5] 

dT 

i^-^ = V{^VT)+jE, (1) 

the Maxwell equation 

and the related equation of the critical state 
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Here u = v{T) is the heat capacity, k = k{T) is the thermal conductivity coefficient, jc is 
the critical current density, and jV is the current density in the resistive state. 

The above system is essentially nonlinear because the right-hand part of Eq.(l) contains a 
term describing the Joule heat evolution in the region of a resistive phase. Such a set (l)-(3) 
of nonlinear parabolic differential equations in partial derivatives has no exact analytical 
solution. 

Let us consider a planar semi-infinite sample (x > 0) placed in an external magnetic field 

-> dH 

H = (0, 0, He) growing at a constant rate = const. According to Maxwel's equation (2), 

there is a vortex electric field E = {0,Eg,0) in the sample, directed parallel to the current 
density j: E || j; where Hg is the amplitude of the external magnetic field and E^ is the 
amplitude of the external electric field. 

For the automodel propagating waves of the type ^{x,t) = x — vt the initial differential 
Eqs. (l)-(3) acquire the following form 
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E = -H. 
c 



(6) 



A set of boundary conditions corresponding to Eqs. (4)- (6) is as follows: 

dT 

T(e ^ +oo) = To, — (e ^ -oo) = 0, 

a? 

E{^ ^ +oo) = , E{^ ^ -oo) = Ee , 



(7) 
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Let us use the Bean-London model of the critical state for the dependence jc{T, H); = 
[6]. According to this model the dependence jc{T) can be represented as jc{T) = jo[l ~ 
a(T — To)]: where jo is the equilibtium current density, a is the constant parameter and 
To is the temperature of the cooling medium. A characteristic field dependence of jr{E) in 
the region of sufficiently strong electric fields {E > Ef) can be aproximated by a piecewise 

linear function jV ~ o'fE, where ctj is the effictive conductivity. In the field of small field 

E 

strengths {E < Ef) we assume a relationship jr{E) ^ Ji In — ,to be valid with ji being 

Eq 

a characteristic local current density scatter (related to a pinning force inhomogeneity) on 
the order of ji ^ 0,01jc and Eq being a constant. This relationship between jr and E is 
due to a thermoactivated early experiments (see,e.g.[l]). The thermoactivated flux creep is 
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principal mechanism of the energy dissipation during the magnetic flux penetration deep 
into the sample. This factor must affect the character of the nonlinear thermomagnetic wave 
propagation in the system studied. 

Excluding the variables T'(^) and H{^), using Eqs.(4) and (6) , and taking into account the 
boundary conditions (7), we obtain an equation describing the electric field E{^) distribution: 



dz"^ 
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The dimensionless variables used in Eq. (8) are defined as follows: 



E.. 



K 



where L is the depth of the magnetic fiux penetration into the sample, is the differential 
conductivity. The corresponding equation of state is obtained using a relationship 



= x{E) 
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where 17= ^ . 

1i^^T2^E^ 

As seen from the plots of XiE) versus E (see figure), there exists a single point of inter- 
section of the curves y = QE^ and y = 1 + (ji/jo) ln{E/Eo), which corresponds to a single 
stable equilibrium state. The stability of this state is determined by the sign of derivative 

d'^E/dz'^ in the vicinity of the equilibrium point. The wave velocity ve can be determined 
dE 

from stationary (— — = 0) Eq.(8) with an allowance for the boundary conditions [7]: 
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Now we can use Eq. (6) and readily derive an equation describing the field distribution 
for a nonlinear H - wave. The wave velocity vh of this H- wave is given by the formula 
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As seen, the wave velocity vh exponentially increases with the field amplitude Hf.. For a 
suffuciently small amplitude He < Ha, where 

-2tjoE,^'/^ 



H„ 



(12) 



the wave velocity is negigible small, which corresponds to the case of the thermoactivated 
flux creep. For H^ = Ha, the wave propagates at a finite constant velocity. In this case, 
the maximum heating of a superconductor in the region immiadiately in front of the wave 
is described by the relationship 

T-Tn H^ 



To 



(13) 
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In conluding, we must note that taking into account a nonlinear dependency of the current 
density j on the electric field strength E would not qualitatively change the main results, 
since the character of the equilibrium state on the phase plane remains esentially the same 
(see [2]). 
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